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ABSTRACT 

We examine high signal to noise A'MM-A''ew;tonEuropean Photon Imaging Camera (EPIC) and 
Reflection Grating Spectrometer (RGS) observations to determine the physical characteristics of the 
gas in the cool core and outskirts of the nearby rich cluster A3112. The XMM-Newton Extended 
Source Analysis Software data reduction and background modeling methods were used to analyze the 
XMM-NewtonEPlC data. From the EPIC data we find that the iron and silicon abundance gradients 
show significant increase towards the center of the cluster while the oxygen abundance profile is 
centrally peaked but has a shallower distribution than that of iron. The X-ray mass modeling is based 
on the temperature and deprojected density distributions of the intra-cluster medium determined from 
EPIC observations. The total mass of A3112 obeys the M — T scaling relations found using XMM- 
Newtonand Chandra observations of massive clusters at r5oo. The gas mass fraction fgas = 0.1491;q q32 
at r5oo, is consistent with the seven-year WMAP results. The comparisons of line fluxes and flux limits 
on the Fe XVII and Fe XVIIl lines obtained from high resolution RGS spectra indicate that there is 
no spectral evidence for cooler gas associated with the cluster with temperature below 1.0 keV in the 
central < 38" 52 kpc) region of A3112. High resolution RGS spectra also yield an upper limit to 
the turbulent motions in compact core of A3112 (206 km s~^). We find that the energy contribution 
of turbulence to total energy is less than 6 per cent. This upper limit is consistent with the amount 
of energy contribution measured in recent high resolution simulations of relaxed galaxy clusters. 

Subject headings: X-rays: galaxies: clusters 



1. INTRODUCTION 

Galaxy clusters represent the largest scales of orga- 
nized matter in the Universe. The intra-cluster medium 
(ICM) traces the cluster gravitational potential created 
primarily by dark matter. It emits X-rays due to the 
highly ionized gas which has been heated by infall from 
the intergalactic medium and enriched by early super- 
novae. Although hydrostatic models produce tempera- 
ture and density profiles in good agreement with obser- 
vations, in the most centrally condensed cases the core 
emission rates are high enough to radiate a substantial 
fracti on of the gas therm al energy over the cluster life- 
time ( Fabian et al.|[T994 ). Metal ions {e.g. iron, silicon, 
and oxygen) recombme as the temperature falls, resulting 
in large increases in the specific cooling rate and produc- 
ing an effective run-away to a nearly neutral gas below 
10'* K and are associated with burst of star formation 



( McDonald et al.| 2011). Assuming pressure equilibrium, 
a continual infall of gas into the core accompanies the 
temperature drop, which over the cluster lifetime, would 
be responsible for a large fraction of the cD galaxy mass. 
This general picture, dubbed as a "cooling flow", was 
also inferred in smaller organized structures (groups) and 
even large elliptical galaxies. Observations with the Ein- 
stein Observatory showed almost a third of all X-ray clus- 
ter s show ed significant apparent cooling fiows (White et 
allp97|. 



nder standing what happened to the cool gas in cool- 
ing flow clusters was one of the prime problems in extra- 
Electronic address: [ebulbuWcfa. harvard. ed"u| 



galactic astrophysics until analysis of the XMM-Newton 
Reflection Grating Spectrometer (RGS) spectra showed 
that they do not, in fact, cool d own to the runaway tem- 



perature (Peterson et al. 2003|. This result shifted at- 



tention to more difficult question: "What terminates the 
coo ling flow?" A numb er of models have been proposed; 
see [ Peterson fc Fabian] ( [200 6 ) for a general review of un- 



derstanding ot cluster core observations and models 

From short early observation s of putative cooling flow 
clusters, Peterson et al. (2003) inferred that there was 
a 'fioor' to the temperature of the cooHng gas at ^ 1/3 
the temperature of the gas, Tmax, just outside the 'cool- 
ing radius' where the cooling time was on the order of 
0.2 or less of the Hubble time. The analysis of RGS 
spectra found less than 10 per cent of the cooling rates 
expected from the observed X-ray surface brightness pro- 
files. Recent analyse s of hi gh-quality spectra of coohng 
flows (Sanders ct al. 2011; Takahashi ct al. 2009; SandersI 



et al.pOOS, , Werner ct al. 2006a. b, Peterson et al.^200: _ 
have shown that the distribution of emission measure 
versus temperature does not have a universal form as ex- 
pected in constant pressure cooling flows and that there 
is indeed gas in some systems over a much wider temper- 
ature range. The wide temperature range in the Centau- 
rus cluster shows that there exists gas with a cooling time 
of less than 10^ years ( Sanders et al.|20d8 |. On the other 
hand i n HCG 62 there is a very narrow range of tempera- 
tures (Sanders et al.'2010). The physical reason for these 



differences is as yet unknown and may provide the miss- 
ing piece to the cooling flow puzzle. The abundances of 
the elements in some of these high signal-to-noise spec- 



2 



TABLE 1 
XMM-NewtonObservations of A3112 





Instrument 


Obs Id 


R.A. 


Decl 


Exposure Time 


Clean Time 


Obs Date 








(J2000) 


(J2000) 


(ks) 


(ks) 




XMM-Newton 


RGS, EPIC 


0603050101 


03 17 57.4 


- 44 14 14.7 


118.9 


102.2 


2009 July 21 


XMM-Newton 


RGS, EPIC 


0603050201 


03 17 57.4 


- 44 14 12.8 


80.8 


80.5 


2009 Aug 08 



tra are also unusual, for example in Centaurus the ratio 
of oxygen (O), magnesium (Mg) and nickel (Ni) to iron 
(Fe) (with respect to solar value) grows monotonically, 
from 0/Fe~l/3 solar to Ni/Fe^2.5 solar value. These 
extreme non-solar abundance patterns may hol d clues to 
the nature of star formation in these systems (O'Dea et 
"ar]|2008 ) and the origin of the energy that prevents run- 
away cooling. 



dra and Suzaku observations of A3 11 2 ( Nevalaincn et al. 
20031 i Bonameiite et al]|2007j [Lehto et ar.||2010^ . |Bona7 



mente et ai. (|2007[) argue that the excess .X.-ray emis 




Fig. 1.— XMM-Newton image of A3112. The white arrow 
shows the dispersion direction for RGS observations. The 5' scale 
corresponds to 410.25 kpc at the cluster's redshift. 



A3112, at a redshift of 0.075, is a rich X-ray cluster, 
wit h a powerful rad io source, PKS 0316-444, in the cen- 
ter (i Takizawa et"aL|,2003k AS CA (Edge et al 1992 ) and 
ROSAI' ( [Allen & l''abian|1997[ ) observations showed that 
A3112 has a strong apparent cooling flow, with an in- 
ferred mass deposition rate of M ~ 415 Mqyi^^. As 
implied by ASCA results, the bright central cD galaxy 
ESQ 248-G006 is inte racting with the surrounding ICM 
( Takizawa et al.|[2003 |. They also detected disturbances 
in the 1(JM around two radio "lobelike" regions near the 
center of A3112 (r ~ 10") from Chandra observations, in- 
dicating tha t the X-ray gas and ra dio lobes are interact- 
ing. Finally, Takizawa et al. ( 2003 1 also found a temper- 
ature decrease and abundance increase toward the center 
of the cluster and placed an upper limit to the contribu- 
tion of X-ray emission from a cooling flow of 10% of the 
total emission. 

In addition to these results, a soft X-ray excess emis- 
sion has been reported from XMH-Newtonand Chan- 



sion is of a non-thermal nature and can be explained 
with the presence of a population of relativistic electrons 
with ~7% of the cluster's gas pressure. They also noted 
that this soft excess emission is equally well fit by a non- 
thermal power-law model or by a thermal model of <0.62 
keV g as temperature and zero metal abundance. |Lehto| 



et al. (2010) showed that the excess emission is ^^10% 
of the cluster hot gas emission at 0.4—1.7 keV assum- 
ing that the emission above photon energies of 2 keV 
originates purely from the ICM based on Suzaku and 
XMM-Newton observations. 

In this paper we examine the results of the deep (182 
ksec) XMM-A^ewtonEuropean Photon Imaging Camera 
(EPIC) and the RGS observations of A3112. This pa- 
per is organized as follows. §2 describes the XMM- 
NewtonFiPlC and RGS data processing and background 
modeling methods. In S|3] we determine Fe, Si, S, and O 
abundance profiles and masses of A3112 out to rsoo based 
on the EPIC spectra products. In ^ we present the re- 
sults of RGS observations. We place an upper limit to 
the turbulent gas motions and the mass deposition rate 
as a function of temperature. We discuss the implications 
of our results and provide our conclusions in Sj5] Errors 
correspond to the 90% confidence level throughout the 
paper unless otherwise noted. In the mass calculations, 
we assume the cosmological parameters h = 0.73, flM= 
0.27 and f^A = 0.73. 

2. OBSERVATIONS AND DATA ANALYSIS 

2.1. XMM EPIG Observations 

A3112 was observed with XMM-Newtonior 119 ks in 
2009 July and 80 ks in 2009 August (see Table [l]). 
The EPIC data processing and background modeling 
were carried out with the XMM-NewtonExtended Source 



Analysis Software {XMM-ES AS) and m ethods ( [Kuntz 
fc Snowden||2008[ [Snowden et"al1|2008[ ). We processed 
the MOS data with the XMM-NewtonScience Analysis 
System (SAS) version 10.0.0 and PN data with version 
11.0.0. The basic filtering and calibration were applied 
to event files using XMM-ESAStools epchain, emchain, 
mos-filter, and pn-fiUer. The tasks pn-filter and mos- 
filter call the S'^S'task espfilt to filter the data for ob- 
vious soft proton flares (SP), using the filtered events 
files created by emchain and epchain. The task espfilt 
creates two light curves for the MOS: one in the field of 
view (FOV) and one from the unexposed corners of the 
instrument in a high energy band from 2.5 keV to 12 
keV. Then the high-energy count rate histogram, which 
should have a roughly Gaussian profile (possibly with 
a high count rate tail) was created from the FOV light 
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Fig. 2. — Light curve and light curve histograms of MOSl obtained from observation 0603050201. The top panel shows a Gaussian 
shaped curve fit to the histogram where there is flaring characteristics due to soft proton contamination. The blue vertical lines indicate the 
range for the Gaussian fit, while the red vertical lines show the upper and lower bounds for filtering the data. The green curve shows the 
Gaussian fit to the data. The middle and bottom panels show the intervals excluded due to soft proton contamination (in black) in the 
FOV and corner light curves. 



curve. The task fits a Gaussian to the peak of the dis- 
tribution and removes flares which exceed a threshold 
count rate set at ±1.5cr level, using the best-fit Gaussian 
parameters. This was used to create good-time-interval 
file. In Figure [2] we show the time periods removed (in 
black) in our analysis due to SP flaring. The exposure 
time after filtering is shown in Table [T] 

Spectra were extracted in a series of contiguous, con- 
centric annuli centered at the X-ray centroid after the 
point sources were removed using the automated point 
source detection XMM-ESAStas\i cheese-hands. The 
products were examined visually to correct for the central 
point sources that were missed by the task cheese-hands 
due to the separation limit we used in the task. The 
center of the image was determined using xmmselect in 
detector coordinates for each observation. We used the 
XMM-ESASta&ks mos-spectra and mos-hack to create 
the spectra in the 0.3—10.0 keV energy band and model 
the particle background and instrument response for the 
MOS detectors, pn-spectra and pn-back were used to cre- 
ate PA^ spectra and background files in the 0.4—10.0 keV 
energy band. The resulting MOS and PA^ spectra of the 
innermost 30" region are shown in Figure ^ 

Proper treatment of the background, along with cross 
talk effects, are the keys to XMM-Newton observations 
of extended sources. The proper methods used in our 
XMM-Newton MOS and PN data analysis are described 
below. 




1 2 

Energy (keV) 

Fig. 3. — An example fit MOSl (turquoise), MOS2 (red), and 
PN (blue) spectra with residual plot from the innermost 30" region 
of A3112 obtained from the observation 0603050201. The diagonal 
lines represent the contamination due t o res idual SP background. 
The fit was performed as described in §2.1| with the model back- 
ground subtracted. The best fit temperature and abundance are 
3.30 ± 0.02 keV and 1.08 ± 0.08 respectively. 



2.1.1. Instrumental Background 

The instrumental background consists of two compo- 
nents: instrumental fluorescent lines and SP background. 
The instrumental fluorescent lines Al Ka (E ~ 1.49 keV) 
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and Si Ka (E - 1.75 keV) in the MOS and the Al Ka 
line (E ^1.49 keV) from the PA^ spectra were modeled 
with two Gaussian components. 

Although the hght curve fihering was applied during 
reduction process as described in §2.1[ some SP residuals 
may still be present in the data. To minimize the effect 
of SP residuals we model the residual SP contamination 
with a single power law model using diagonal spectral re- 
distribution matrices obtained from XMM-ESASCmrent 
Calibration Files database ^. The power-law index was 
left free to vary between 0.2 and 1.3 with a free normal- 
ization factor for the innermost spectra. The normaliza- 
tions of the power-law model of the other annuli were tied 
to the normalization of power-law of the innermost re- 
gion according to their area. The power-law background 
models for MOSl (black), M0S2 (red), and PN (green) 
data are also shown in Figure [s] for the innermost 30" 
region. 

2.1.2. Cosmic X-Ray Background 

The cosmic X-ray background (CXB) is spatially vari- 
able over the sky and has multiple components. The 
CXB was modeled by two thermal components for the 
emission from the Local Hot Bubble (LHB) and the 
Galactic halo, and a non-thermal component for the un- 
resolved cosmic sources. We use the ROSAT All— Sky 
Survey (RASS) spectrum to model the CXB using the 
background tool at the High Energy Science Archive Re- 
search (HEASARC) Web site. The RASS spectrum was 
extracted from an annulus surrounding the cluster center 
between 1° and 2°, with the assumption that the annulus 
spectrum reasonably represents the CXB in the direction 
of the cluster. We simultaneously modeled the emission 
from the LHB (or heliosphere) with a cool unabsorbed 
single temperature thermal component (E ^ 0.1 keV) 
and the Galactic hotter halo and intergalactic medium 
with an absorbed thermal component (E ^ 0.2 keV). 

Another CXB component arises from the emission 
from unresolved point sources. The contamination due 
to unresolved cosmological sources is modeled using an 



absorbed power law com 
a ^-^1.46 (Snowden et al. 



jonen t with a power law index, 



2008) 



TABLE 2 

O VII and OVIII emission surface brightness 

Instrument Line Energy Line Surface Brightness 
(keV) (LU) 



MOS 


OVII 


0.56 


3.92 ±1.05 


MOS 


OVIII 


0.65 


1.75 ± 0.56 


MOS-I-PN 


OVII 


0.56 


-I +0.61 
^- ' ' -0.72 


MOS-I-PN 


OVIII 


0.65 


1 66 

1-DD -0.61 


2.1.3. 


Solar 


Wind Charge 


Exchange 



The cosmic soft X-ray background below ^1 keV is 
dominated by line emission from oxygen ions, such as the 
0.56 keV O VII (the O VII triplet from n = 2 ^ 1) and 
0.65 keV OVIII (the O VIII Ly a transition). The origin 




0.0001 



Energy (keV) 

Fig. 4.— 60 MOSl, M0S2 and PJV spectra obtained from ob- 
servations list ed in Table[T]were fit simultaneously using a method 
described in i ]2.l| with a HASS spectrum. The total reduced 
obtained from this fit is 1.28 for 14650 degrees of freedom. We plot 
five MOSl and a RASS spectra surrounding the A3112 centeroid 
obtained from observation 0603050201 as an example. 

of these lines could be the hot gas in the LHB, charge ex- 
change between Oxygen ions in the solar wind (SWCX) 
and geocoronal or interplanetary materi al, or a combi- 



nation of the these two processes, (e.g., Snowden et al 



Smith et al. 2007 Fujimoto et al. 2007|). The poten 



i2004j 

tial S WCX emission was modeled by two Gaussian com 
ponents for lines of OVII at E - 0.56 keV and OVIII at 
E ^ 0.65 keV with widths set to zero. To distinguish be- 
tween the LHB X-ray emission and the SWCX emission, 
we calculate line fluxes in units of photons cm~^ sT^ sr~^ 
(hereafter LU, for "line units"; see Table [2]). While rela- 
tively large, these values are not outside the range seen in 
other observations; see, e.g., Koutroumpa et al. (2011). 
Without high-resolution non-dispersive spectroscopy de- 
termining the true origin of this emission is not possible. 

2.1.4. Cross Talk 

Due to the finite point spread function (PSF) of 
the XMM-Newton'EPlC detectors, significant numbers of 
photons originating from one region may be detected in 
another region if the region annuli are comparable to the 
size of the PSF. This is an especially important source 
of bias in temperature m easurements of strong cool core 
clusters, such as A3112 
null widths of 0.5' and 



2008). The an- 



il Snowden et al. 

we used in this work require 
the use of a crosstalk correction as an additional model 
component in the spectral analysis. The cross-talk effect 
leads the measured temperature of the innermost annu- 
lus to be higher and of the outer annulus to be lower than 
the true value. We used the SASiask arfgen to create 
cross-talk ancillary region file (arf) to correct for pho- 
tons originating in one annulus and detected in another 
annulus. 



2.1.5. 



EPIC Spectra Fitting 

The fitting w as done i n the spectral fitting package 

^XSPEC 12.6.0 ( |Arnaud|[l996| . The 0.3-10 keV en- 

^ http : // xmm.vilspa.esa.es /external /xmm_sw-cal/calih/ccf.shtml\ eT:gy band was used for MOS fits whereas 0.4—10 keV 
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Fig. 5. — 1st order and 2nd order RGS spectra and background fit with a single temperature thermal model. The best fit temperature 
and abundances are 3.82 ± 0.21 keV and 0.78 ± 0.02 respectively. The of the source spectra is 4149.5 for 2732 de grees-of-freedom and 
background fit is 2167.1 for 1332 degrees-of-freedom. 



TABLE 3 

Elemental Abundances Obtained from EPIC Observations. Upper 
limits are 90 per cent. 



Region 


Fe 


Si 


O 


S 


(arcsec) 


(Solar) 


(Solar) 


(Solar) 


(Solar) 


- 30 


l.llitO.Ol 


0.94±0.07 


0.82±0.04 


0.87±0.09 


30 - 60 


0.54±0.02 


0.12±0.11 


0.53±0.08 


< 0.17 


60 - 120 


0.41±0.02 


< 0.11 


0.37±0.12 




120 - 180 


0.32±0.03 




0.36±0.15 




180 - 240 


0.34±0.04 




0.57±0.22 




240 - 300 


0.35±0.05 




0.36±0.23 




300 - 420 


0.24±0.04 




0.19±0.21 





band was used for PN fits. The Galactic absorption col- 
umn density was left free. The spectral fitting was done 
simultaneously for a total of 61 spectra using the ex- 
tended C— statistics. The diffuse cluster X-ray emission 
was fit with an absorbed single- temperature colli sional 
equilibrium plasma model (apec) ( Smith et al.|2001| ) and 
AtomDB v2.0.1 (A. Foster et al. 2011, in preparation). 
The final value was obtained from the best fit deter- 
mined by C-statistics. Since C-statistic does not provide 
a direct estimate of the goodness-of-fit and computing 
the null hypothesis probability distribution of simulta- 
neously fit 61 spectra requires a large amount of compu- 
tational power over a long period of time we assessed the 
goodness of the fit using the value we calculated from 
the final fit obtained from C-statistics. We found that to- 
tal reduced was of 1.28 for 14650 degrees of freedom 
(dof) from this simultaneous fit. The temperature and 



abundance parameters were constrained to be the same 
for MOSl, M0S2 and PN detectors but allowed to vary 
for each region. A bundances were given with respect to 
the solar values in Anders &: Grevesse ( 1989 ) . The nor- 
malizations of the single temperature thermal model for 
the cluster emission were linked for the MOS and PN 
observations. The temperature for the CXB model was 
linked for different regions and instruments. The OVII 
and OVIII line parameters representing the SWCX were 
linked for MOS and PN data. Figure |4] shows the five 
MOSl spectra and a RASS spectrum obtained from the 
final fit as an example. 

2.2. RGS Observations 

The high resolution grating spectrometer (RGS) pro- 
vides high-resolution X-ray spectroscopy in the wave- 
length range from 5 to 38 A. We examined two RGS 
observations (0603050101 and 0603050201, see Table [l]) 
of A3112 using the version 10.0 of the Science Analysis 
Subsystem {SAS ) software. The SAS task rgsproc was 
used to create the filtered event file and source spectrum. 
Time periods where the background count-rate exceeds a 
threshold count-rate value 0.2 counts s~^ were removed. 
The light curves of background events created from the 
filtered event files showed no sign of SP flaring. The 
exposure times after the filtering process are shown in 
Tabled 

RGS spectra were extracted by making selections both 
in the spatial (dispersion/cross-dispersion) and energy 
(dispersion/PI) planes. The default 90 per cent of the 
PSF of a point source selection region in the cross- 
dispersion direction was used for the extraction of source 
spectra including 90 per cent of the peak of the pulse- 
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Fig. 6. — Radial profiles of temperature and abundance for A3112 obtained from the single temperature thermal model (apec) fit. The 
temperature and abundance measurements were obtained from the MOSl, MOS2, and PN data are shown in purple. The MOSl and 
MOS2 results are indicated in red. The blue data points are the measurements reported in |Snowden et ar] | |2008] |. 



height distribution. We generated the RGS background 
and response matrix files using the S'^S'tool rgshkgmodel 
and rgsrmfgen. The RGSl spectra obtained from two ob- 
servations (0603050101 and 06030501201) fisted in Table 
[T] were co-added using the task rgscombine. The same 
process was applied for RGS2 observations. The com- 
bined RGS 1st order and 2nd order spectra were then fit 
in the 7 — 28A and 7 — 16A wavelength range respec- 
tively (see Figure [5]). 

Due to the slitless nature of RGS spectrometers the 
observed line emission from an extended source is broad- 
ened according to the spatial extent of the source. The 
observed broadening in wavelength (AA) is 



AA = ^A.A 



(1) 



where A9 is the source extent in arc-minutes and m is 
the diffraction order. ^ The slitless RGS spectrometers 
do not allow one to distinguish a change in X-ray photon 
wavelength and a change in position along the disper- 
sion direction. The most straight-forward way to correct 
for the line smearing due to source extend is to use SAS 
routine rgsrmfsmooth prior to spectral fitting. The rou- 
tine modifies the response matrix by convolving it with 
a spectral model based on the angular structure function 
computed from a high-resolution MOS or Chandra im- 
ages in the RGS dispersion direction. This approach is 
based on a simple approximation that the X-ray spec- 
trum is uniform a nd the spatial distribut ions of all the 
lines are identical ( Rasmussen et alT||2001 ). 



3. EPIC RESULTS 

The temperature and abundance profiles, shown in 
Figure [6j obtained by fitting 61 MOS and PN spectra 
are obtained from both observations (see Table [I]) with a 
single temperature thermal model (vapec). The oxygen 
(O), silicon (Si) and iron (Fe) and sulfur (S) metallicities 
were allowed to be free in the fit with the other elemen- 
tal abundances tied to Fe, except He which was fixed to 
solar. 

^ \http : / /xmm.esac.esa.int /external /xmm.user_support/\ 
documentation / uhb I rgsspecres.html 
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Fig. 7. — Radial distributions of Iron (Fe), Silicon (Si), Sulfur 
(S) and Oxygen (O) abundances with 90 per cent confidence limits 
obtained from the XMM-Newton MOS and PN observations. The 
best fit values are obtained by fitting MOS and PN spectra with 
an absorbed vapec model. 

We compare our temperatu re and abundance mea sure- 
ments with those reported by Snowden et al. (2008) (red 
data points in Figure [61) . Our temperature ana abun- 



dance profiles obtained from combined MOS and PN fits 
are consistent with those measured by [Snowde n et al.| 
(1^08). The difi'erences in the outskirts is mainly due to 
the exposure time di f ferenc e between this work (200 ksec) 
and [Snowden et al.| (2008[ ) (20 ksec). M05 temperature 
measurements are systematically higher than MOS+PN 
combined temperatures shown in the left panel of Fig- 
ure [6[ A systematic difference of up to 7% in the 
wide band temperatures of other clust ers between MOS 
and M OS+PN detectors was noted by |Nevalainen et al.| 
(20101. The abundance measurements obtained from 
MUS observations are in good agreement wi th those ob- 
tai ned fro m MOS"-!- PA^ observations and the [Snowden et| 
( |2008|r esults. 
Figure |7[ shows the radial distribution of Fe, Si, S and 
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O elemental abundances obtained from the MOS and PN 
spectra. The fit was performed with an absorbed single 
temperature thermal model (vapec). The Fe, Si and O el- 
emental abundances were allowed to vary independently 
while carbon (C), sulfur (S) and argon (Ar) were coupled 
with Fe abundance as we were unable to obtain good 
constraints for the abundance of these elements. We find 
that the Fe abundance shows a significant increase to- 
wards the center of the cluster while the O abundance 
profile is centrally peaked but has a shallower distribu- 
tion than that of Fe. Si was detected only in innermost 
60". Nevertheless, there is also an indication of an in- 
crease in the Si abundance in the cluster core. Type la 
supernovae (SNe la) produce significant amounts of Fe 
and Si while type II supernovae (SNe II) produce large 
quantities of O but very little Fe and Ni. The centrally 
peaked trend of Fe is qualitatively consistent with the 
idea that the contribution from SNe la relative to that 



from SNe II increases towards t he cluster center ( Kaastra 
et aLj2001[|Tamura et al.|2004[ ). Fe is still being added to 



the f CM specifically in the core by the SNe la, and O is a 
well-mixed product of the SNe II that explode soon after 
the first star bursts. The centrally peaked distribution 
of O may be due to ram-pressure stri pping of gas from 
in-falling galaxies ( de Plaa et al.|[2007| . The comparison 
of metal abundances obtained from FIPIC and RGS ob- 
servations with current SNe models will be addressed in 
the upcoming paper (G. E. Bulbul & R. Smith (2012, in 
preparation). 

3.1. Radial Profiles of Cluster Masses 

In order to examine the cluster properties as a function 
of radius, we obtained the projected density from the 
normalizations of the best-fitting thermal plasma model 
to MOS data we used in XSPEC. The normalization is 
related to the emission measure as 



Norm ■ 



10- 



-14 



47r(Z)A(l + z))2 



nenndV, 



(2) 



where rie ~ 1.2 for a fully-ionized cosmic plasma and 
Da is the angular diameter distance. Instead of correct- 
ing for projection by combining a series of models, we 
use a model-independent approach to subtract the pro- 
jected spectra from each s uccessive annulus to produce a 
set of deprojected spectra (|Fab ian et al.|198l" Kriss et al. 
19831. This direct 'onion peeling' deprojection method 



assumes a spherical geometry and that the density is con- 
stant within each annulus. The projected normalizations 
of the single temperature thermal model (APEC) were 
deprojected using the 'onion peeling' method. These de- 
projected normalizations were then used to generate the 
deprojected number density profile using Equation [2] 

We describe the density distribution of the hot plasma 
in galaxy cl usters using the model proposed by |Bulbul| 



et al. (2010). The deprojected density profile was repre- 



sented by 



1 {l + r/rsf-^ -1 
{P-2) r/r,{l + r/r,f-^ 



cool 



Ur) (3) 



where rigo is the normalization factor for the density pro- 
file and n is the polytropic index. The core taper function 
Tcooh used for cool core clusters such as A3112 is 



'^cool 



{r) = 



a + {r/rcooi)'^ 
1 + {r/r^ooiV 



(4) 



where rcooi is the cooling radius, a and 7 determine the 
st rength and shape of the cooling in the cluster center. 



ture profile, which is 



Bulbul et al. (20101 model also provides the tempera- 



T{r) = To 



1 (l + r/rs)'^-^ -1 
(/3-2) r/r,(l + r/r,)'3-2 



oi{r). (5) 



where Tq is the normalization factor for the temperature, 
Ts and (/3 — 1) are the scaling radius and the slope of the 
total matter distribution. 

In order to fit the temperature profile obtained from 
MOS observations, we projected the three-dimensional 
temperature model shown in Equation [5] along the line- 
of-sight. This projection approach is based on a weight- 
ing method that correctly predicts a good approximation 
of the spectroscopic temperature for a mixture of differ- 
ent te mperature and abundance plasma fMazzotta et al. 
2004) rViklilinin„2006( ) 
ture T ' 



The spectroscopic-like tempera- 



is, 



Jy w r-3/2 dV 



(6) 



where w is the weighting factor w — n'^ j^-a ^^^^ ^ 
^ 0.75 for the XMM-N ewtonand Chandra observations 
( Mazzotta et al.|[2004[ ). The spectroscopic-like temper- 
ature obtained using the three dimensional temperature 
profile (Equation [s]) was then compared to the tempera- 
ture profile obtained from MOS observations. The result 
of this fit is shown in left panel of Figure [8| 

The fitting was performed using the Monte Carlo 
Markov Chain algorithm, with Metropolis-Hastings sam- 
pling, to determine posterior distributions for the best fit 
model parameters. The best-fit parameters and 90 per 
cent confidence intervals of the model are shown in Table 
[4} We obtained a = 13.3 for 12 dof for the combined 
temperature and density profile fits. 

The gas mass Mgas is computed by integrating the gas 
density profile within the volume, 



Mgas{r) = 47r^e rup / n^ir) r^ dr 



(7) 



where /ig and nip are the mean molecular weight and the 
proton mass. The total mass is 



M(r) = 



Anpir^ 



1 



1/(1 ~ /3) - r/r. 



(/3-2) V/9-1 (l + r/r,)'9-i 



Tcoolir). 

where fi is the mean molecular weight and nip is the pro- 
ton mass. The normalization factor for the total matter 
density is = (ToHn + l)(/3 - l))/iATrGpmprs^). The 
gas mass fraction is 



fgas 



Mr, 



(9) 



We measure gas mass, total mass and gas mass fraction 
at two physically meaningful overdensity radii, 7'25oo and 
rsoQ. The overdensity radius ta is defined as the radius 
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Fig. 8. — Deprojected electron number density and temperature profiles for A3112. Left panel: the number density profile where the 
black points are derived from the normalization of the single temperature thermal fit from XMM-Newton MOS spectra, the red line and 
the shaded area shows the best fit model and the 90 per cent confidence limit obtained from the model fit. Right panel: the black data 
points are the temperature measurements obtained from MOS spectra. The red line shows the best fit spectroscopic-like temperature 
profile obtained from the projection of Bul bul et al.| |2010 ) model and the shaded region shows the 90% confidence level obtained from the 
model fit. The combined of the fit is 13.3 for 12 degrees of freedom. 

TAB LE 4 

Best-fit Parameters of the \Bulbul et al. \ 1(201 Of Model 



(10^2 cm-^) 


(arcsec) 


(kpc) 


'^cool 

(arcsec) 


r-cool a 7 
(kpc) 


To 
(kcV) 


1 17+0.66 1 oq+0.29 Q cn+1-68 
^•^'-0.32 ^■°'^-0.27 '^•"'-'-0.78 


64 Q+48.8 


00 7+66.7 
°°-'-35.8 




115.5i^|;^ 0.03 0.73t!^;J^ 





TABLE 5 

Gas Mass and Total Mass Measurements at Overdensity Radii i?2500 '^'^'^ R500 



TA Mgas Mtot 

(arcsec) (kpc) (IO^Mq) (lO^^j^^^) 



fgas 



Reference 



2500 331.0i«i 452.5l26;4 1.68«;}? 1.53l«;2| O.lOQig;^^^ 



2500 



459 



1.7 ± 0.30 



This work 



Nulsen et ah (20101 



500 708.8 
500 



+43.7 

36.5 



8.9 

1025 



+59.7 
68.2 



4.49 



+0.43 
0.45 



J-UU_o 44 

9 Q+2.4 



0.149; 



-0.036 
-0.032 



This work 



Nulsen et al. I 2010 1 



within which the average matter density of the cluster is 
A times the critical density of the universe at the cluster 
redshift: 



(10) 



where Pcrit is the critical density of the Universe at the 
cluster's redshift. The measurements of gas mass, to- 
tal mass and gas mass fraction at two overdensity radii, 
'''2500 and rsoo, are shown in Table [5] We find that 
the total mass measurements obtained from our mea- 
surements are in ag reement with the previous studies 
(Nulsen et al. (2010); see Table [S] for comparison). We 



also estimate the total mass using the Mtot — T scal- 
ing relations derived from the XMM- iVewtonobservations 



of ga laxy clusters (Arnaud et al. 2005 Mantz et al. 
20101. We find that the average temp erature within 
0.15 V boo < r < rsoo is ~4.6 keV. The [Arnaud et al.| 
(|2005| Mtot — T scaling relation produces a total mass 



oOI^ot = (3.24 ± 0.12) X IO^^Mq for this temperature, 
which is consistent with our result at the l-cr level. The 

T scaling relations produce a 



Mantz et al. 



J2010D Mtot 

total mass oi" M^= (3.43 ± 1.82) x IO^^Mq and are also 
consistent with our result at the l-cr level. 

Figure|9 shows the gas mass fraction profile. Our result 
agrees with the Wilkinson Microwave Anisotropy Probe 
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Fig. 9. — Gas mass fraction profile of A3112. The dashed hori- 
zontal line indicates the WMAP measurement of the cosmic baryon 
budget Qt/^m = 0.167 l |Komatsu et al.120111 . 

(WMAP) m easured cosmic baryo n fraction flh/ilM = 
0.167at rsoo ( [Komatsu et al.||2011[ ). 



4. RGS SPECTRAL RESULTS 

The default 90 per cent of the pulse-height (P.I.) distri- 
bution was used for the extraction of source spectra. The 
diffuse X-ray emission from the first-order and second- 
order RGS observations were then fit with an absorbed 
single-temperature coUisional equilibrium plasma model 
(apec) and AtomDB 2.0.1 (A. Foster et al. 2011, in prepa- 
ration). The best-fit temperature, abundance and the 
goodness of the fit is given in Table [6) Since C-statistic 
does not provide a direct estimate of the goodness-of-fit 
and we asses the goodness of the fit using the value we 
calculated from the best-fit obtained from C-statistics. 
The final values reported in Table [g] were obtained 
from the best fit determined by C-statistics. 

4.1. An Upper Limit to the Source Extent 

The two different region selection filtering was used 
for the extraction of source spectra (the default 90 per 
cent and 95 per cent). The best-fit temperature, abun- 
dance and the goodness of the fit for these two different 
aperture are given in Table [6) In order to determine the 
source extent observed by the RGS we determined the 
line width of the O VIII Lya line without applying the 
spectral broadening convolution model to the response 
files for two different pulse height extraction mask. The 
best fit model to data for O VIII Lya observed line en- 
ergy (0.607 keV, 20.44A) and width (2.26x10-^ keV, 
0.043A) is shown in Figure [lo| The observed line emis- 
sion has been broadened by a combination of the thermal 
broadening of the ICM, turbulent motions within the ob- 
served region and the spatial extent of the source along 
the dispersion direction. Assuming that the O VIII line 
is intrinsically narrow and only broadened by the spa- 
tial extent of the source we placed an upper limit to 
the source extent observed by RGS to be ~ 38" (~ 52 
kpc) using Equation [l] This result shows that the RGS 
is only sensitive to emission from the cluster core, and 
retains useful resolution in these observations. The tem- 
perature and metal abundances measured by the RGS 
spectra are consistent with the central bin temperature 
obtained from MOS and PN observations (see Figure [6]) . 
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Fig. 10. — Close up view of oxygen VIII Lya line obtained from 
high resolution XMM-Newton RGS observations of A3112. The 
best fit Gaussian line is shown in blue. The background spectrum 
is shown in red with the best fit model (dashed red line). This fit 
produces total 64.3 for 50 degrees of freedom for source spec- 
trum and 89.8 for 43 degrees of freedom for background spectrum. 

4.2. Elemental Abundances in the Core 

We fit the RGS first order and second order spectra 
with a single temperature thermal model {vapec). The 
O, Ne, Mg, Si and Fe abundances were allowed to vary 
while all other metal abundances were tied to Fe abun- 
dance. The normalizations were allowed to vary between 
first second and second order spectra. The best fit tem- 
perature and elemental abundances obtained from this 
fit are shown in Table [7] with the goodness of the fit. 

We find that the temperature of the gas in the cluster 
core obtained from iZCS* spectra agrees with the temper- 
ature of the central bin (0 — 30") of EPIC observations. 
The abundance measurements of Fe, O and Si in the in- 
nermost 38" region obtained from RGS observations are 
consistent with the abundance in the central bin of EPIC 
observations at 1-a level (see Figure [7]). 

4.3. Multi-temperature Fit 

In order to test if there is a multi-component gas in 
the core as found in other clusters, e.g ., Centaurus (|Taka- 



hashi et al.|2009[ [Sanders et al |2008| , we fit i^GS" spectra 



with two vapec models each with a variable temperature, 
normalization. The Ne, O, Mg and Fe metallicities were 
allowed to vary independently while the Si, Ca, Al, S , Ar 
and Ni metallicites were linked to Fe and tied between 
two different vapec components. 

We also extended the two-temperature model to six 
temperature components with fixed temperatures ranges 
4.0, 2.0, 1.0, 0.5, 0.25 and 0.0808 keV to fit the RGS 
spectra. We coupled all metal abundances for each model 
to reduce the number of free parameters but allowed each 
normalization to vary. Table [7[ shows the best-fitting 
parameters for these six vapec models. The temperature 
weighted emission measure for each component is plotted 
as a function of temperature in Figure [TT] Figure [Tl] 
indicates the apparent lack of any gas below 1 keV in 
the RGS results. A3112 is particularly different from 
the previous results re ported for Centaurus , HGC 62, 
A262 and A3581 (Sanders et al.[[2008[ [2010| ) since RGS 
observations of these clusters show a multiphase gas at 
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TABLE 6 

RGS-fit Results for Two Different Region Selections 



P.I. Width* 


Source Extent 


kT 


Abundance 


Red. 


(per cent) 


// 


(kcV) 




(d.o.f.) 


90 


37.9 


3.82 ± 0.21 


0.78 ± 0.02 


1.36 (3287) 


50 


30.6 


3.52 ± 0.37 


0.86 ± 0.17 


1.60 (2468) 



*The inclusion region of the pulse- liciglit mask. 



TABLE 7 

Temperature and Abundance Best-fit Parameters of RGS Data Extracted From the 90 per cent PSF aperture. 



Parameter 


vapec 


2 X vapec 


vapec + vmcflow 
kT^^ji free 


6 X vapec 


kT (keV) 


3.89 ± 0.26 


4.75 ± 0.25 


4.78 ± 0.33 


4.0, 2.0, 1.0, 0.5, 


kT„i„ (keV) 




1.01 ± 0.03 


0.52 ±0.07 


0.25, 0.0808 


O (Aq) 


0.74 ± 0.14 


0.90 ± 0.13 


0.72 ± 0.12 


0.75 ± 0.03 


Ne (Aq) 


2.10 ± 0.40 


2.52 ± 0.35 


2.36 ± 0.30 


2.21 ± 0.21 


Mg (Aq) 


1.10 ± 0.26 


1.45 ± 0.30 


1.46 ± 0.37 


1.70 ± 0.20 


Fe {Aq) 


0.90 ± 0.12 


1.33 ± 0.11 


1.57 ± 0.30 


1.28 ± 0.024 


Mass deposition rate 






43.5 ± 3.5 




Norm (x 10"*) 


127.7 ± 0.17 


119.5 ± 2.3 


94.5 ± 6.4 


See Figure [ll] 






2.76 ± 0.55 






(d.o.f.) of source spectra 


4149.5 (2732) 


4067.0 (2730) 


4122.0 (2731) 


4039.6 (2727) 


(d.o.f.) of background spectra 


2167.1 (1332) 


2382.4 (1332) 


2410.9 (1332) 


2370.9 (1332) 



lower temperatures (^0.5 keV). 
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Fig. 11. — Emission measure as a function of temperature ob- 
tained from RGS observations. RGS spectra was fit with six vapec 
models with fixed temperatures of 4.0, 2.0, 1.0, 0.5, 0.25, and 0.0808 
keV. The best fit abundances and goodness of the fit are shown in 
Table [7l 



4.4. Fitting the Cooling Flow 

We also determine the cooling rate of the gas as a func- 
tion of temperature using the latest AtomDB v2.0.1. The 
XSPEC model vmcflow models gas cooling between two 
temperatures and gives the normalization as a mass de- 
position rate. Here the mass deposition rate is an upper 



limit to the gas which cools below X-ray temperatures 
determined from the X-ray spectra. We constructed a 
model composed of a vapec model and a cooling flow 
model components (vmcflow) to fit the i?G'5' spectra. We 
left the kTmin parameter free in order to determine the 
minimum temperature that gas is cooling to whereas the 
maximum temperature was linked to the temperature of 
the vapec model. The abundances of each element was 
allowed to vary while the redshift was fixed to the value 
reported in Table [7] The mass deposition rate we ob- 
tained from this fit is 43.5 ± 3.5 Mq yr~^ . We also 
found that there is no evidence for the emission below a 
cut-off at ~1 keV. The lack of emission from cool gas in- 
dicates that this emission property is general for clusters 

with a cD galaxy. 

From Einst ein ( [Edge et al.|1992| and ROSAT ( [Allen fc 
Fabian 1997) imaging data with little spectral informa- 
tion a large total mass dep osition rate has been reported 



for A3112 of 415 M© yr'^ O'Dea et al. (2008) reported 



that the total mass deposition rate was much less than 
previously claimed with 10^5 Mq yr~^ being consistent 
with the Chandra and XMM-Newton X-ray spectra. The 
level of cooling we report in this work is also consistent 
with the star formation r ate of 4.2 Mp yr~ ^ measured 
from the Spitzer IR data ( [O'Dea et al.||2008 |. 



4.5. An Upper Limit on Turbulence 

The turbulent gas motion in the ICM can be gener ated 
by frequent mergers ([Ascasibar fc Markevitchl2006[) , the 
motion of galaxies ( [Conro y & Ostriker 2008), SJN explo- 
sions and, AGN outbursts ( jMcNamara fc IMulscn 2007). 
The line widths of heavy elements are dominated by ther- 
mal broadening. Very low (~50— 150 km s"^) turbulent 
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TABLE 8 

Measured Value of the total velocity broadening. If line smearing 
effect due to the spatial extended of the source is ignored this 
velocity corresponds to turbulent velocity. 



TABLE 9 

90 Per cent Upper Limits to Turbulent Velocity Broadening 
(^turb.\\)- rgsrmfsmooth was used for removing the effect of the 
line broadening due to the extended nature of the source. 



Line 


(km s^-"-) 


Energy Fraction 
(%) 


Line 


(km s~^) 


Energy Fraction 

(%) 


OVIII line-only 






OVIII line-only 


417 


24 


bvapec (all lines) 




54±2 


bvapec (all lines) 


206 


6 



motions can affect the magnetic field orientation in the 
ICM and increase the effectiv e level of thermal conduc- 
tion ( Ruszkowski fc O"h||2010 ). The dissipation of turbu- 
lent gas motions has been proposed as one of the mech- 
anisms that heats the ICM in the core of clusters and 



prevents t he gas from extrem e cooling (Rebusco et al. 
[2008 ; Ruszkowski fc Oh|[20lol ). 

The earliest direct upper l imit on turbulence in the X- 
ray waveband was placed by Sanders et al. (2010) exam- 
ining the XMM-NewtonRGS spectra, of A1835. I'he limit 
on the total broadening, including the turbulent compo- 
nent, of 274 km s~^ was reported at the 90 per cent 
level. It was also reported that the ratio of turbulent to 
thermal ener gy density in th e core of A18 35 is less than 



13 per cent. Sanders et al. (2010 2011) placed upper 



limits on the velocity broadenmg examming the XMM- 
NewtonRGS observations of 62 galaxy clusters, galaxy 
groups, and elliptical galaxies. After the addition of a 
50 km systematic uncertainty it was reported that 
five of these sources to have less than ~20 per cent of 
the thermal energy density is related to turbulence, with 
velocity upper limits <220 km s~^. 

We examine the deep RGS observations of A3112 and 
can place an upper limit on velocity broadening of the 
heavy metals in the ICM. The total line width due to 
thermal motions of ions and turbulence is, 



'turb ' 



'the 



(11) 



the line of sight component of line width due to turbu- 
lence is 



1^0 



Vturb,\\ 



(12) 



and of the line of sight component of line width due to 
thermal motion of ion i is 



rrij 



(13) 



Here Vq is the rest frequency of the line of interest, Te 
is the electron temperat ure, nij is the mass o f the ion i, 
and c the speed of hght ( Rebusco et al.|[2008 ). Then the 
line of sight velocity of turbulence is 



'^turb.l 



1/2 



(14) 



The energy density contribution due to turbulence can 
be calculated via, 



f-the 



fllTlp 



^turb,\\ 



(15) 



where fi is the mean molecular weight, Wt„rb || is the line 



of sight velocity, and nip is the proton mass (Werner et 
,al 2009. ). 

in order to test an extreme case we first assume that 
the line broadening due to source extent is negligible and 
lines are broadened purely by the thermal motions of 
ions. The best fit total broadening obtained from the 
RGS spectra was used in Equation [14] to calculate the 
velocity of turbulence. The energy contribution of tur- 
bulence to the total energy was calculated using Equation 

110 



15 



We found that the velocity of broadening is 804 
km s^^ with an energy contribution of 52^2*' P^r cent by 
examining only OVIII line region between 18.5 A and 22 
A wavelength range. We also fit RGS first order and sec- 
ond order spectra with bvapec model which utilizes the 
AtomdDB 2.0.1 atomic data base and includes a velocity 
fit paramet er that accounts for line broadening ( Sanders 
The entire wavelength band (7 — 28A) we 



et al. 



2010) 



used for this analysis gives a measurement of 8341 I^q^ km 

s-i (see Table [8|. 

We then correct for the line broadening due to the 
source extent as described in {2.2 by convolving response 
matrices with the MOS image. This fit produces a total 
of 4061.4 for 2650 dof for the source spectra and 2866.1 
for 1537 dof for the background spectra. We find a 90 
per cent upper limit to turbulence velocity of 206 km 
s~^ with an energy density contribution of 6 per cent for 
the entire i?G'>S' wavelength band (7 — 28A). Using only 
the OVIII line (18.5 A and 22 A) wavelength range we 
obtain a 90 per cent upper limit of 417 km s~^ with a 
energy contribution of less than 24 per cent (see Table 
^and Figure 10). Total we obtained from this fit is 
64.3 for 50 doTior the source spectrum and 89.8 for 43 
dof for the background spectrum. 

We repeat the analysis of modifying the response ma- 
trices convolving with a Ghandra image obtained from 
the offset AGIS-I observations to correct for the line 
smearing due to the spatial structure of the source. We 
obtain a total of 4003.1 for 2650 dof from the source 
spectra from this fit. Background spectra produce a to- 
tal of 2840.0 for 1537 dof. This analysis produces 
a 90 per cent upper limit to the velocity broadening 
of 225 km s~^ for the entire RGS wavelength band (7 
— 28A). The vignetting and off-axis effects of Chandra 
AGIS-I image could partially be responsible for the dif- 
ference between upper limits to the velocity broadening. 
For self-consistency in the analysis we report upper lim- 
its obtained using the MOS image as described in j ]2.1[ 
We also note that rgsrmfsmooth routine which is used to 
correct for broadening due to the spatial structure of the 
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TABLE 10 

90 per cent Upper Limits to Line Fluxes and Metal Abundance obtained from RGS Spectra 



Line 


Transition 


^rest 


Flux 


Abundance 






(A) 


( 10^" photons cm^^ s^^) 


< 


Fe XVII 


3d -s> 2p 


15.014 


1.00 


4.98x10"" 


Fe XVII 


3s -s> 2p 


16.780 


0.84 


6.89x10"'* 


Fe XVII 


3s -s> 2p 


17.096 


6.97 


4.25x10-3 


Fe XVIII 


3d 2p 


14.208 


5.45 


5.58x10-3 


Fe XVIII 


3s -5- 2p 


16.071 


2.40 


4.39x10-3 



source makes an approximation that the source spectrum 
does not change as a function of spatial position. 

Non-thermal pressure support from bulk motions in 
clusters cores plays an important role in the thermal 
properties of galaxy clusters. In particular, the Sunyaev 
Zel'dovich (SZ) effect observables for cluster surveys, 
such as the Atacama Cosmology Telescope, the South 
Pole Telescope and Planck, CARMA/SZA are affected 
by these processes (|Nagai et al. 2007[ Shaw et al 



Battagha et al.|20Hp . Several earlier simulations s. 



2010 



tiowed 

that the kinetic pressure from bulk motions contributes a 
small ( ~ 5% —15%) but still significant amount of energy 
within cluster core scales in rel axed clusters, which are 
consisten t with out upper limit ( Dolag et al.|2005 Lau et 
al.|2009 )■ More recent works, performing high-resolution 
simulations with the AMR code ENZO simulations of 
galaxy clusters have shown that the energy contribution 
of turbulent motions, Eturb/ Etherm, is <5 per cent in the 
compact cores of relaxe d clusters, are consistent with our 
upper limit |Tap ichino fc Niemeyer|2008 Lau et al.|2009 
Burns et al.|[5niOf |Vazza et al.||201l| : 



4.6 



Absence of Low Temperature Component and An 
Upper Limit on Metal Abundance 

The soft excess was found earlier in A3112 in the 



XMM - Newt o n, Chandra and Suzaku d ata by Nevalainen 
et al.| ( |2003| ), [Bonamente et al.| ( |2007[ ), and |Lehto etliL 
( |2010 ). Interpretation of the soft excess emission has 
been under debate. One possible explanation is that this 
excess emission originates from a thermal w arm gas as- 
socia t ed with the cluster (kT ~ 0.1—1 keV) (Lieu et al. 
1996 Bonamente et al. 20071, which indicates the pres- 
ence of an additional ph ase m the intergalactic medium. 
Bonamente et al. (2007) have placed an upper limit on 
metal abundance of < 0.01 by fitting the soft excess emis- 
sion detected from the XMM-MOS observations with an 
additional single temperature thermal model. 

We, therefore, searched for the signatures of < 1 keV 
gas in the RGS spectra in order to test the hypothesis of 
the thermal origin of the soft excess detected from this 
cluster and to determine a new upper limit on the metal 
abundance by assuming that the soft excess follows the 
cluster's surface brightness. Since the Fe XVII and Fe 
XVIII lines are sensitive indicators of a plasma with a 
temperature of < 1 keV, we determined the Fe XVII line 
flux limits, are shown in Table [TOl To put a limit on any 
emission from Fe L shell lines, we added a 6 function to 
the model at the energy where the Fe XVII and Fe XVIII 
lines would be expected. The 90 per cent upper limits 
to the Fe XVII and Fe XVIII fluxes are shown in Table 
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Fig. 12. — 90 per cent upper Limits to the Fe XVII and Fe 
XVIII Line Fluxes obtained from the RGS observations A3112 are 
shown in green, red and blue straight lines . The curves show the 
expected line fluxes for [Bonamente et al.| p007J temperature of 
0.62 keV plasma, shown m vertical black Ime, with normalization 
of 1.05x10-^ obtained from AtomDB 2.0.1 database (A. Foster et 
al. 2011, in preparation). 



the limits we report are for compact emission from Fe 
XVII, XVIII lines. The limits on Fe XVII, XVIII lines 
from extended emission are presumably larger. 

These upper limits were then compared with the line 
fluxes estimated using AtomDB v2.0.1. The fluxes 
were estimated assuming an optically-thin low-density 
plasma. Best-fit values for the MOS observations from 
Bonamente et al. ( 2007 ) were used to convert from emis- 
sivities to line fluxes: the geometric normalization factor 
of 1.05x10-^ and column density of 2.6 x 10^" cm"^. 
The emissivity of the Fe XVII and Fe XVIII ions peaks 
in the 0.5—0.7 keV temperature range, with predicted 
fluxes 1 to 3 orders of magnitude higher than the ob- 
served 90 per cent upper limits obtained from RGS spec- 
tra at the pl a sma t emperature 0.62 keV given by [Bona^ 
mente et al. (2007), as shown in Figure 10 The upper 
limits determined from Fe XVII and Fe XVill line fluxes 
are shown in Table [TO] This implies that the product 
of the abundance and normalization (AxN) must be 100 
times less. These abundance and normalization products 
we obtained from RGS spectra allow us to place a new 
upper limit on the metal abundance. These upper lim- 
its we obtained from i^GS" spectra indicate that 0.62 keV 
plasma associated with the cluster is either very tenuous, 
to the point of non-existence, or has a very low metallic- 
ity less than 4.98 x 10~* in the central < 38" region of 
the cluster. 
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5. CONCLUSION 

In this work we studied the X-ray emission of a nearby 
rich cluster A3112 from deep XMM-Newton EPIC and 
RGS observations. The XMM-Newton EPIC data pro- 
cessing and background mode hng were carried out with 



the most up-to-date methods ( Kuntz & Snowden 2008 
Snowden et ah 2008). The teniperature and abunidance 



profiles obtained from our combined MOS and PA'' fits 
are co nsistent with those measured by [Snowden et aL| 
(2008). We find that MOS temperature measurements 



are systematicahy higher than M06'-|-i-'iV combined tem- 
peratures. This systematic discrepancy in the wide band 
temperatures of up to 7 per cent was also observed by 
Nevalainen et al. (2010). The abundance measurements 
obtained from the MUS observations are in good agree- 
ment with th ose obtained from the MOS+PN observa- 
tions and the Snowden et al. (2008| results. From the 



EPIC fits we find that the iron abundance gradient shows 
a significant increase towards the center of the cluster 
while the O abundance profile is also centrally peaked 
but has a shallower distribution than that of Fe. The 
centrally peaked trend of Fe is qualitatively consistent 
with the idea that the dominant contribution from SNc 
la is towards the cluster center, whil e the distribution of 
SNe II remains more or less uniform ( Kaastra et al.|200T 
Tamura et al.||2004|. This result may also indicate that 



iron is still being added to the ICM specifically in the 
core by the s upernovae SNe la. 



We use the Bulbul et al. (2010) model to examine the 
cluster properties as a function of radius. The projected 
density and temperature are obtained from the MOS 
spectra. We then deproject the density profile obtained 
from XSPEC normalizations and project the three di- 
mensional temperature model to fit the MO S data. We 
demonstrated that the Bulbul et al. (2010) model ac- 
curately describes the number density and temperature 
profiles from the obtained from the XMM-NewtonX-Ta,y 
data. At r^oo, the total mass M500 = (3^00l^J|)x 



lO^^M p) , is in agreement with the results of Nulsen et al. 
(2010). The total mass also obeys the empirical M — T 



scal ing relations found i n general massive ga laxy clus 



ters 



Arnaud et al] ( |2005[ ); [Mantz et al.| ( |2010[ ). The gas 

mass fraction is fgas = 0.1491q;o32 ^500i and is consis- 
tent with the seven-year WMAP results (Komatsu et al. 



We also searched for the signatures of < IkeV gas. 
The comparison of flux limits on the Fe XVII and Fe 
XVIII lines in the RGS sp ectra with the At omDB line 
fluxes is obtained using the Bonamente et al. (j2007j) val- 



ues for the plasma temperature, column density, abun- 
dance, and the geometric normalization factor. The up- 
per limits on the Fe XVII and Fe XVIII lines we obtained 
from RGS observations allow us to place an upper limits 
on metal abundance in the central 38" region of A3112. 
These upper limits indicate that either there is no spec- 
tral evidence for cooler gas associated with the cluster 
with temperature below 1.0 keV in the central 38" (< 52 
kpc) central region of A3112 or either it has a very low 
metallicity (< 4.98 x lO"'^). 

We also determined the emission measure distribution 
as a function of temperature, detecting or limiting the 
amount of cool gas in the cluster core. Alternatively, 
the emission should have a cut-off at ~1 keV. The lack 
of emission from cool gas has been found also in similar 
systems , such as A1835 (|P eterson~et 



al 



200T|), Sersic 



159-(3 3 ( [Kaastra et al.||200ip , and A1795 ( jTamura et ah 
2001 1, indicating that this emission property is general 
for clusters with a cD galaxy. 

We also placed an upper limits to turbulent gas mo- 
tions directly using high-resolution i?C?S' X-ray spectra of 
the compact core of A3112. The upper limit of 206 km 
s^^ to turbulent motions with less that 6 per cent con- 
tribution to the total energy that we found in the core of 
the relaxed cluster A3112 is compatible with the values 
of ~6 per cent of pressure in gas motion s, calculated us 



ing high-resolution Eulerian simulat ions ( Lau et al.]"2009_ 
Burns et al.|2010 Vazza et al.|2011 ). These upper limits 



were obtained based on an assumption that the source 
spectrum does not change as a function of spatial posi- 
tion. The low limit we measured on turbulent motions 
from the core of a relaxed galaxy cluster A3112 is en- 
couraging for the cosmological studies using SZ cluster 
surveys. In the near future, the launch of Astro-H will 
provide an X-ray calorimeter with high energy resolution 
of ^4.5 eV. This detector will produce measurements of 
lines widths with high accuracy revolutionizing the mea- 
surement of turbulent velocities in the ICM. 
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